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Abstract 

A method for determining the leading quantum contributions to the effective 
action for both zero and finite temperatures is presented. While it is described in 
the context of a scalar field theory, it can be straight-forwardly extended to include 
fermions. An extrapolation procedure which can significantly enhance the compu- 
tational efficiency is introduced. This formalism is used to investigate quantum 
corrections to the nucleation rate in first-order phase transitions. 
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1. Introduction 



The evaluation of quantum corrections to classical solutions is an important problem 
which pervades much of modern theoretical physics. However, while effective potentials have 
been studied extensively, methods for determining the effective action are less well-developed. 
Moreover, the actual evaluation of such effects for realistic systems has often been hampered 
by their general intractability to analytical solution and the lack of efficient computational 
methods. t 1-5 ' 

In the effective potential approximation to the effective action, quantum fluctuations 
are integrated out about a constant classical field — but this is not expected to be adequate 
because the classical field is generally an inhomogeneous configuration. The derivative ex- 
pansion^ improves on this by accounting for spatially varying background fields; its leading 
term is the effective potential. The expansion is a perturbative approximation which ex- 
tracts the dominant contribution of short-distance quantum effects on long-distance physics. 
When it converges, it provides an efficient means for performing calculations. However, 
when it diverges, one must often resort to brute- force techniques which entail an explicit, 
computationally-intensive evaluation. Furthermore, the derivative expansion fails whenever 
the potential V is non-convex (V" < 0) in some region of space, which includes an important 
class of perturbatively calculated potentials.' 7 ' It is clear that a general method, which is 
also applicable to such cases, is needed. 

In this paper, a method for calculating the quantum effects arising from the effective 
action is presented. The next section contains the general formalism for evaluating the 
effective action. Section 3 discusses the exact formulation of the computational method as 
well as extrapolation techniques which improve its convergence properties. This formalism 
is applied to the analysis of phase transitions in Section 4J 7 ' 8 1 



2. General Formalism 

Consider a scalar field theory with the Lagrangian density 

£(0) = ^^-^), (2.1) 

where V is the tree- level potential which has a (classical) vacuum at <p v . Since this paper 
considers 3+1 dimensions exclusively, renormalizability constrains V to have no polynomials 
in cj) of higher power than a quartic. The classical field <p is determined by the equation 

d li P$=-V'(4>). (2.2) 
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The contribution of one-loop quantum effects to the effective action can be written as 



Tr£ (1) = -Trln 
2 



8 2 + V"(4>) 



d 2 + ^ 2 



where /x 2 = V"(4> v ) and the trace runs over space-time coordinates. t Part of this trace can 
be evaluated as TrC™ = Tr' / (t\C^\t)dt, where Tr' runs over the remaining spatial degrees 
of freedom. Since this paper deals only with time-independent fields, specializing to this 
case means that states in the energy basis \co) are eigenstates of the operator in So 
inserting a complete set of such states and performing a partial integration yields 



Tr£ (1) = -i / Tr' 



-uj 2 - v 2 + v"(4>) -uj 2 - v 2 



2 — dt. 
2tt 



Observe the non-locality of this expression; this generic feature of loop corrections makes 
exact analytical treatments difficult. The remaining trace can be conveniently performed 
over the eigenstates of the operators in C^: if ip® and ^ - are chosen such that 



[-V 2 + ^? = (a/J)VJ, 

[-v 2 +n0)]^ = K) 2 ^, 

where the subscript j indexes the eigenstates, then 

^ (1) = ~/EK-- w ?)*- 

Hence the one-loop effective action can be written as 



(2.3) 
(2.4) 



S cff (0) = J [£($)- £(&)] d*x + f 



„°<A 



dt. (2.5) 



The bare sum in TrX^ is divergent; it is regulated in eq. (2.5) by a momentum cut-off A, 
and a counterterm £ ct ((f>, A) has been added to render it finite. 

For time-independent fields 0, it is more convenient to focus on the energy E of the 
system which is related to S eS through 



S, 



cff 



(0) = ~E(4>) J 



dt. 



^The trace excludes possible negative and zero modes of the operator d 2 + V"{<j>). When such modes arise, 
they must be explicitly removed and treated differently. I 7 ] 
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Then 

E{4>) = E cl (4>) + E^((f), A) + E ct ((p, A), (2.6) 
where E cl is the energy of the classical field configuration 

EM = ~J [m - d 3 x, (2.7) 

E<U is the one-loop contribution 

£ {1) (0,A) = i ]T(^-^), (2.8) 

w o<A 

and £J ct is the energy due to the one-loop counterterms 

E ct (^A) = -J CM,A)d 3 x. (2.9) 

At finite temperature T, the free energy F replaces E:^ 

F cl = E ch AF 1 = E^\ F ct = E ct ; (2.10) 

for bosons, <f>(x, r) is periodic in Euclidean time r with period T _1 , and there is an additional 
contribution due to one- loop effects:' 7 ' 

AF T = T^( T — Jpf ) (2.11) 

Observe that no additional counterterms need to be added to F ct because finite temperatures 
do not change the short-distance behaviour of the theory. 

In the following section, we describe the method developed to evaluate the the quantum 
corrections, AF 1 and AF T , formally given by eq. (2.8) and (2.11). While for some special 
situations, the w ■ can be obtained analytically, this is unfortunately not possible for a general 
potential V"(4>). Instead the eigenvalues must be found numerically, then for AF 1: the bare 
sum ^2 u ,°<a( uj i ~ ^ s computed explicitly, and finally the counterterm subtracted; for 
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AF T , the sum in eq. (2.11) must be performed term- by-term. To attain reasonable accuracy 
this subtraction has to be done at a large cut-off A (to achieve convergence) when both 
the bare sum and the counterterm (which individually diverge as a function of the cut-off) 
are numerically very large. Since the final result is much smaller, each term has to be 
determined very precisely, resulting in a heavy computational burden. Furthermore, the 
straight-forward approach of evaluating the free energy by a "brute-force" term-by-term 
summation of the expressions in eq. (2.8) and (2.11) until convergence is reached is also 
computationally inefficient. 
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3. Method of Computation 



To circumvent the above- described problem of having to compute both the regulated 
bare sum and its counterterm to very high numerical accuracy, the three-dimensional problem 
is first decomposed into channels of definite angular momentum. Then for each channel, 
the divergent part of the bare sum is analytically removed through subtraction with the 
corresponding divergence in the counterterm, leaving a much smaller finite piece. Since the 
contribution of higher partial-wave channels decrease rapidly, this procedure overcomes the 
problem. 

An improved computational method is then presented. It is based in part on the obser- 
vation that the higher-energy modes in the spectrum of eq. (2.4) are less perturbed by the 
potential V"(<f>) due to the non- uniform background field </> than the lower-energy ones. This 
allows us to formulate an approximation method which accounts for the contribution of the 
high-energy modes accurately (where the accuracy of the approximation increases with the 
energy) so that only some of the lower-energy modes need to be treated exactly' 2 ] 

Exact Formulation 

The difference in the eigenenergies uj^ and uj® of the unbound states (uj > /i) can 
be characterized by the phase shift between the (asymptotic forms) of the corresponding 
continuum state eigenfunctions ipj and ip®, as was first shown in one dimension.' 1 ] Since 
the phase shift is generally a well-behaved, smoothly-varying function of the energy, it is 
relatively easy to calculate. Hence it is convenient to express the free energy in terms of 
this quantity. To determine the phase shift, we consider eq. (2.4) which determines the 
fluctuations about the classical field configuration. 

Since most classical solutions <fi exhibit spherical symmetry (0 = 4>(r)), we will restrict 
our analysis to such systems. Then the solution to eq. (2.4) can be separated into radial and 
angular parts by choosing an eigenfunction of the form 




(3.1) 



where the radial wavefunction is determined by 




(3.2) 



with the boundary condition u n l (0) = 0. The Y lm are the spherical harmonics corresponding 
to a state with total angular momentum I and ^-component m. 

The corresponding equation for u° nl where V"{4>) is replaced by fi 2 , 




(3.3) 
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has an exact analytical solution: 



u° nl (r) = V2k nrJl (k n r) (3.4) 

where j l is a spherical Bessel function and k 2 = oo 2 — ^ ■ These solutions have the asymptotic 
form 

In 

u° n i(r)^V2sin(k n r-—),r^oo. (3.5) 

The potentials we consider behave asymptotically as V"(tj)(r)) — > fx 2 when r — > oo (which 
corresponds to those with finite action). For such potentials, the asymptotic behaviour of 
the solution to eq. (3.2) will be 

In 

u n i( r ) -> V2sin{k n r - — + 5,(fcJ), r -> oo. (3.6) 

These equations serve to define the phase shift ^ for each angular momentum channel I. 
Note also that both u n l and l are (2/ + l)-fold degenerate. 

To facilitate the counting of states, it is convenient to discretize the eigenvalue spectrum. 
This procedure can be achieved by enclosing the system in a box of radius L (where L is 
much greater than the range of the potential V") and imposing the boundary condition 

u n j{L) = (3.7) 

which requires that 

Note that such a discretization is implicit in the formal sums in eq. (2.8) and (2.11). The 
values attained by cu° (before discretization) as defined by the energy eigenvalue of eq. (3.3) is 
a continuous spectrum ranging from an energy of fi to infinity. The corresponding spectrum 
for u determined by eq. (3.2) will generally consist of some discrete bound states with 
energies iv? < fx 2 and a continuous spectrum with energies u 2 > fi 2 . The difference in 
structure between the continuum spectra of the two systems manifests in a difference in the 
respective density of states. Hence it is appropriate to express the sum over eigenenergies 
for the states in the continuum as an integral over the density of states: 

POO 

J^wJ = ^(21 + 1) / ojrf(cu)duj (3.9) 

l J u 



5>= E 



poo 

(21 + l)co nl + ^2(21 + 1) / uj ni (uj)duj (3.10) 

i J u 
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where (21 + 1)7^(0;) is the density of states of angular momentum / for the potential V"(4>) 
with an analogous definition for On taking the continuum limit (L — > oo), it follows 
from eq. (3.8) that the densities of states are related to the continuum phase shift through 



7r da; 



(3.11) 



Now if eq. (2.4) has N bound states, then since eq. (2.4) and eq. (2.3) must have the same 
total number of states, 



N + J2(2l + 1) n l (uj)dLU = ^2(21 + 1) I 

l J ii l J n 



n^(u)du. 



(3.12) 



For a finite potential, this implies Nit = 5(fj). 

It is convenient to define the free energy in each angular momentum channel such that 



AF 1 (A) = ^(2/ + 1)AF{(A), 
i 

AF T = J^i 21 + 



(3.13) 
(3.14) 



and to similarly partition the counterterm energy as 



F ct (A) = ^(2/ + l)Ff(A), 



(3.15) 



then from the above equations 



/ ,.,2 ^,,2 Z7F ; -V 



u))duj, 



and 



AF? = T £ In (L 



e-^/ T / 7r /., e^/ 7, - 1 



da;. 



(3.16) 



(3.17) 



The Appendix discusses the renormalization of the Lagrangian given by eq. (2.1). It is 
shown there that the contribution of the counterterms to the energy are of the general form 



J g(u 2 -p 2 -^) J h{x)S 



X 



d 3 p duj 
(2tt) 3 2^ 



/CO J 
T*(9h)^, (3.18) 



where g is a power of the propagator, h is a function of 4> and its derivatives, and Tr' is 
a trace over the spatial variables. The partial wave decomposition of these contributions 



6 



is achieved by taking the trace with respect to the eigenstates of eq. (2.3) denoted here by 
\nlm): 



Tr'(gh)=J2 (nlm\g{u 2 + V 2 - fi 2 )\n'l'm')(n'l'm'\h{r)\nlm) 



nlm n'l'm' 
i$>+l) 



/ 2 2 2\ 

(u -p - n ) 



h(r)\u° pl (r)\ 2 dr 



dp. 



(3.19) 



From the Appendix, the counterterm contribution to the free energy is 



F ct (A) = -- f Tr'U (u,p) [m 2 (r) -/i 2 ] + \a o (u, P ) 2 [m 2 {r) - /i 2 ] 



21 2 



duo 



(3.20) 



where 



A (W,P) = — ~ 2 J—": 

ur — p l — n z + le 



and m 2 (r) = V"{4>{r)). Evaluating the trace using eq. (3.19) yields 



1 



1 



\u 



° pl (r)\ 2 [m 2 (r)-^} 



4tt (p 2 + ^ 2 )V2 y 



(3.21) 



where A = a/A 2 — /i 2 is a three-momentum cut-off. 

This completes the formulation of the method for the exact calculation of the free energy. 
However, as we have remarked above, the convergence of such an exact computation can be 
sufficiently slow so that extrapolation techniques can be useful. Amongst the various such 
procedures, we consider in particular the WKB approximation, which provides an analytic 
expression for the phase shift that is valid at high energies and hence can significantly reduce 
the effort required to evaluate the phase shift integral. ^ 

WKB-Improved Method 

A differential equation of the form 



f(x) = 



(3.22) 



has an approximate WKB solution given by 

exp [i /* k(u, y)dy] 



fwKB( x ) 



(3.23) 
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which is valid when the wavelength is much less than the distance scale over which k varies: 



1 dk 



where k(cu,x) is the local wavenumber 



k(u,x) = \Ju 2 - V"(4>(x)). 



Hence the accuracy of the WKB approximation increases with energy. The phase shift for 
such solutions is given by 



rWKB 



» = f 

J — c 



k(u, x) — lim k(u>, y) 

y— >oo 



dx. 



(3.24) 



Explicitly 



^WKB, 



r J CU 2_ m 2 (r) _ i i±ii dr _ r J UJ 2_ fi 2_ i i^ dri ( 3.25) 

Jatui) ' r Ja n (uj) ' r 



where a and a denote the classical turning points defined by 



- 0. and cu 2 - /j, 2 - 



uj 2 — m 2 (a) 



= 0. 



a" 



On 



Applying this method to eq. (3.2) yields an analytic expression for the energy integral of the 
phase shift: 



Jfi Jo 



u 2 - m 2 (r) - l{ - l \ l \ {k - Q{r))dw 

n(r) * r 

l l±lle(A-n (r))du 



A 
Q (r) 



UJ 2 — jJ? 



(3.26) 



dr 



with 



il(r) - V- 2 (r) + ^-^. 



flo(r) 



and is the unit step-function. Observe that since the high-energy behaviour of the 
phase shift is independent of the angular momentum, the energy of each angular momentum 
channel is logarithmically divergent: 



AF'(A) = -L In [m 2 {r) - ^ 2 ] dr + O(A ) 



(3.27) 



Now the divergent piece in AF( can be analytically combined with the infinite part of F t ct in 
eq. (3.21) to leave only finite terms. Performing this subtraction and taking the limit A — > oo 



gives the final expression for the WKB-improved, temperature-independent renormalized free 
energy: 

AF< = lim [Af?(A) + if(A)] 



A-+00 



= 2 E (""-ri-tol 6l{Uj)dU (3.28) 

- ^ { X/ ( A W kb , r) + ^ (r ) [m 2 (r ) - /x 2 ] + (r ) [m 2 (r ) - 2 ] dr 
In this equation, is the contribution from the WKB phase shift above A WKB , 

X,(r) = m " (r)2 ; m ' (r)2 " ^A efr (r)VAe ff (r) 2 -^(r) 2 

+ ^A c ° ff (r) v /A c V(r)2-m['(r)2 

+ i TO|(r) 2 ln A cff (r) + v /A cff (r)2-m,(r)2 (3-29) 
2 ^tx 



where 



--m,(r) in - , 

m°(r) 2 = p 2 + m,(r) 2 = m 2 (r) + ^±H, 

Ae ff (r) = max (A WKB , Q (r)) , A eff (r) = max (A WKB , fi(r)) , 

and A WKB denotes the energy above which the phase shift is computed by the WKB method. 
The remaining terms in the last integral come from finite parts of the counterterm with 

(r)= r s2 \ ^ s )\ 2 -\ d (330) 

iy } Jo V* 2 + (H 2 

and 2 oo 2 2 

ft ( r ) = [°° S \w\ d s. (3.31) 
4 J (s 2 + (pr) 2 ) 3 / 2 

Equation (3.28) indicates that AF^ ren can now be computed by first summing over the 
bound state energies, then the continuum state contributions can be evaluated by explicitly 
computing the exact phase shift only up to A WKB , beyond which the WKB method provides 
an analytical expression that accounts for contributions at higher energies. Note that while 
the WKB procedure entails an approximation, its accuracy can be made such that the 
difference between the exact and the WKB results is smaller than the desired precision. 
Finally, summation over I yields 

AF 1;ren = + !) AF Un- (3-32) 

i 

Since AF T is not divergent, it can be computed exactly using eq. (3.14) and (3.16), or 
by replacing the exact phase shift S t above a certain energy scale by the approximate WKB 
phase shift <^ WKB given by eq. (3.25). 
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4. Application and Discussion 



In ref.[7] these methods have been used to calculate the free energy of an instanton 
configuration which determines the decay rate in a first order phase transition. The compu- 
tation of AF 1 will be described first. It is found that the accuracy available on conventional 
computers prevents a precise determination of this quantity when it is straight-forwardly 
evaluated as in eq. (2.8) — that is, by doing the bare sum and subtracting the counterterm, 
without a decomposition into partial waves. When AF 1 is computed exactly, by utilizing 
such a decomposition, very high numerical accuracy is still required because for each / the 
bare sum and i^ ct (A) must be evaluated at a large value of the cut-off A. But since both 
quantities diverge as a function of A, we find that convergence with reasonable accuracy is 
still difficult to attain. In contrast, evaluation of AF 1 ren using the WKB-improved method 
consisting of eq. (3.28) and (3.32) converges rapidly for much lower values of the cut-off 
A WKB and typically only the first fifty partial waves need to be summed; the parameters 
required for convergence are very much dependent on the nature V"(<j)(r)) and the values we 
have quoted come from the potentials we have examined. 

The exact computation of AF T can be performed by evaluating eq. (3.14) and (3.17), 
but at high temperatures it is found that several hundred partial waves must be summed 
to attain convergence. When the exact phase shift is replaced by the approximate WKB 
expression at high energies, there is a reduction in the computational burden and the same 
number of angular momentum channels must be summed. The improvement is not marked as 
it was for AF l rcn in part because AF T is not renormalized. The results of these computations 
are summarized in Tables 4 and 5 of ref . [7] . 

In summary, we have elucidated a method for the exact evaluation the effective action 
to one-loop. The WKB extrapolation scheme was devised to reduce the computational 
effort. These methods enable an efficient calculation of the free energy associated with a 
phase transition, as detailed above. However, the applicability of this method is not limited 
to this example. Rather, it can be utilized in a broader variety of problems involving the 
non-perturbative evaluation of observables in a non-uniform background in quantum field 
theory! 9 ' as well as in classical systems^ . It can also be generalized to encompass theories 
with fermions.' 11 ' 
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Appendix: One-loop Renormalization of the Scalar Field Theory 



This Appendix discusses the one-loop and renormalization of the scalar field theory 
described by eq. (2.1). The classical vacuum <f> v satisfies 

V'(4> v ) = and f, 2 = V"(cf> v ) > 0. (1) 

At one-loop the only divergent graphs are those with one and two vertices corresponding to 
quadratic and logarithmic divergences, respectively. 

It is convenient to adopt a renormalization scheme where the counterterms are chosen to 
exactly cancel the divergent graphs as shown in Figs. 1 and 2. These conditions are imposed 
at zero external momenta; this choice has the advantage that the one-loop contribution to 
the effective potential V 1 satisfies 

m) = m) = wj = o, (2) 

so that eq. (1) is unchanged at one-loop. Then the counterterm Lagrangian to be added to 
eq. (2.1) is 

4*= a[v(<i>) - ^]+\n v "(<t>) - ^ 2 ] 2 > ( 3 ) 

where 

If i d A k 
a — - 



2 J P- y u 2 + ie(27r) 4 ' 
and 

^~2J {k 2 -^ + ie) 2 (2vr) 4 ' 

The terms in eq. (3) involving a and (3 renormalize the graphs with one and two external 
vertices, respectively. These divergent integrals can be suitably regularized by imposing a 
momentum cut-off A. 
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Figure Captions 

Figure 1. Renormalization scheme for the divergent one-loop graphs with one vertex. A box with 
a cross denotes a counterterm insertion. 

Figure 2. Renormalization scheme for the divergent graphs at one- loop order with two vertices. 
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